TEMPERATURE DEPENDENCE OF THE CASIMIR EFFECT 



I. Brevik^ 

Department of Energy and Process Engineering, Norwegian University of 
Science and Technology, N-7491 Trondheim, Norway 

J. B. Aarseth^ 

Department of Structural Engineering, Norwegian University of Science and 
Technology, N-7491 Trondheim, Norway 



Abstract 

In view of the increasing accuracy of Casimir experiments, there 
is a need for performing accurate theoretical calculations. Using ac- 
curate experimental data for the permittivities we present, via the 
Lifshitz formula applied to the standard Casimir setup with two par- 
allel plates, accurate theoretical results in case of the metals Au, Cu 
and Al. Both similar and dissimilar cases are considered. Concen- 
trating in particular on the finite temperature effect, we show how 
the Casimir pressure varies with separation for three different tem- 
peratures, T = {1, 300, 350}K. The metal surfaces are taken to be 
perfectly plane. The experimental data for the permittivities are gen- 
erally yielding results that are in good agreement with those calculated 
from the Drude relation with finite relaxation frequency. 

We give the results in tabular form, in order to facilitate the as- 
sessment of the temperature correction which is on the 1% level. We 
emphasize two points: (i) The most promising route for a definite ex- 
perimental verification of the finite temperature correction appears to 
be to concentrate on the case of large separations (optimum around 
2 /im); and (ii) there is no conflict between the present kind of theory 
and the Nernst theorem in thermodynamics. 
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1 Introduction 



Consider the standard Casimir configuration, namely two semi-infinite ho- 
mogeneous media separated by a vacuum gap of width a (recent reviews on 
the Casimir effect are given in [U El El EI)- Assume that the surfaces are 
perfectly plane, and that they are of infinite extent. Denote the left hand 
slab by 1, the intermediate vacuum region by 2, and the right hand slab by 3. 
The two permittivities are ei(a;) and £3(0;). Spatial dispersion is neglected. 
Working in terms of complex frequencies we have uj = iC,. The Lifshitz 
variables s and p, and the Matsubara frequencies Cm, are 

q 27im 



P=--, Cm =—5-, (1) 

where (3 = 1/T is the inverse temperature, and m the Matsubara integer 
(we usually put h = c = 1). The nondimensional frequency y and the 
nondimensional temperature 7 are defined by 

y = qa, q = \f^i + CL, 7 = (2) 

k_L being the transverse wave vector (i.e. parallel to the surfaces). Now 
defining quantities A's related to the TE and TM modes by 

eip + si ^sP + ss 

we can write the Casimir pressure as 



1 ^ /-oo 



m- 



=0 •'"^T 



(4) 



where the prime on the summation sign means that the m = term is 
counted with half weight. 

The above formalism summarizes the exposition given recently in j3] . The 
main purpose of the present paper is the following: 

1. We calculate the Casimir pressure more accurately than we did earlier, 
inserting updated permittivity data for the metals Au, Cu and Al (courtesy 
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of Astrid Lambrecht), and extend also the region of separation to larger val- 
ues. The numerical corrections are small, on the 1% level, but the generally 
improving accuracy of pressure experiments accentuates the need for work- 
ing to this degree of accuracy. As before, we take the plates to be perfectly 
smooth; roughness corrections have to be dealt with separately and are not 
covered here. 

A general property of the Casimir pressure is that from T = onwards 
it decreases with increasing values of T, given a fixed value of a. However 
when the separation is large, a > 3 /xm, the pressure increases with increasing 
T when T becomes high. Our numerical results indicate that the pressure 
increases from around room temperature and upwards. As a general con- 
clusion, a separation of about 2 //m appears to be optimal for testing the 
temperature correction to the Casimir pressure; then the relative correction 
is highest (cf. also figure 5 in The practical problem, of course, is that 
at large separations the Casimir pressure itself is small. 

2. We next emphasize the point that there is no conflict between this 
kind of theory (which is equivalent to adopting the Drude dispersion relation) 
and the Nernst theorem in thermodynamics. According to this theorem, the 
entropy at zero temperature has to be equal to zero. The Nernst theorem 
is satisfied in our case, in spite of the fact that the contribution from the 
TE zero mode (that is, the m = term in the Matsubara sum) is equal 
to zero for finite e as well as for a real metal. The point here for a real 
metal is that the relaxation frequency stays different from zero. The first to 
emphasize this kind of behaviour were Bostrom and Sernelius . We have 
treated these thermodynamical issues ourselves also [313131101, and there 
are several other works expressiong the same opinion ^T] Ej. We 

mention, though, that the opposite view has also been advocated in recent 
papers [13 CEl HZl ■ 



2 Calculation and results 

Let us first recall the Drude dispersion relation 

where Up is the plasma frequency and v is the relaxation frequency. The 
plasma wavelength is Ap = l-ncjiOp. or the three metals mentioned, the 
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corrected data as compared with those given in j3] are 



uop = 9.03 eV, u = 34.5 meV, Ap = 137.4 nm Au, 



ujp = 8.97 eV, 



u = 29.5 meV 



A 



138.3 nm Cu, 



Ur, = 11.5 eV- 



u = 50.6 meV- 



A 



107.9 nm Al. 



(6) 



These corrections are roughly on the 1% level. As before, we calculate the 
Casimir pressure by means of MATLAB, extracting the zero-frequency case 
171 = for a separate analytical treatment. Since e becomes very large in the 
zero frequency limit for metals, we can express the m = contribution as 



where the polylog function with arguments (3,1) is involved. 



Calculated values of the Casimir pressure for Al, Cu and Al are shown in 
tables 1-6, both for the similar and the dissimilar cases. As in we took the 
temperature T = 1 K to represent the case T = case with good accuracy. 
In the calculations our tolerance for the integrals was 10~^^, whereas the 
tolerance in the sum was 10"*^. At T = 1 K the necessary number of terms 
was quite large, especially at small separations (for instance, about 25700 at 
a = 0.16 fim). 

It is seen that the room-temperature pressure is always weaker than the 
zero-temperature pressure. Thus for Au-Au plates at separation a = 0.5 /im, 
the pressure is lowered from 16.56 mPa to 15.49 mPa, or by 6.5%. The 
reduction becomes much enhanced at larger separations; thus at a = 2 yum 
the pressure is lowered from 7.549 x 10~^ mPa to 5.550 x 10~^ mPa, or 
by 26.5%. It thus seems advantageous to work with high separations, if 
technically possible. The differences between Au-Au and Cu-Cu pressures 
are generally small, whereas the pressures for Al-Al are larger, as we might 
expect from the dispersive data in 

As for the cases where dissimilar metals are involved, the Au-Cu data, 
table 4, are quite similar to those given in tables 1 and 2. The data for Au-Al 
and Cu-Al, tables 5 and 6, show somewhat larger pressures. 
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1 



polylog(3, 1) 



0.1502571129. (8) 
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In the tables we show also the pressures when T = 350 K, since one 
may expect that the pressure difference between 300 K and 350 K will be 
soon measurable. Again considering Au-Au at a = 0.5 fim, we see that the 
pressure is lowered from 15.49 mPa to 15.30 mPa, or by 1.2%, when T is 
increased from 300 K to 350 K. If a = 2 //m, the corresponding pressure 
decrease is 3.7%. 

A striking property is that for the larger separations, the pressure in- 
creases with increasing values of T. This turns out numerically when a 
becomes larger than about 2.8 /xm. 

Finally, the following point should be noted: we did not have to use the 
Drude relation in any of our calculations, for any finite frequency. All the 
frequencies that were needed, were lying within the region of Lambrecht's 
data. We needed the Drude relation explicitly at only one place, namely in 
the evaluation of the zero frequency term, m = 0. 

3 On Nernst's theorem 

It is important to ensure that the present formalism does not come into con- 
flict with basic thermodynamics. In the present case this means in particular 
that the Casimir entropy per unit area at zero temperature. 



{F being the free energy), has to be zero at T = 0. This is Nernst's theorem. 
Let us make some brief remarks on this topic, in view of its current interest. 

• There exist no measurements of the permittivities at very low frequen- 
cies. What is at our disposal, is a series of measurements of room-temperature 
complex frequencies s{u!) = s'{u!) + e"{Lj), where the data on s"{u!) permit 
us to calculate the real quantities £{iC) via the Kramers- Kronig relation. 
The permittivity data received from Lambrecht cover the frequency region 
1.5 X lO"*""^ rad/s to 1.5 x 10^^ rad/s. Based upon these data, the relaxation 
frequency u in the Drude relation is determined. For low frequencies we 
have to describe the permittivity analytically, with use of the Drude rela- 
tion, down to C = 0. It is not very important, however, to know the value 
oi u a,t T — very accurately; the important point is that u at T — stays 
finite. In practice, this condition is always fulfilled because of scattering from 
impurities. Then, it is easy to show that the zero frequency TE mode does 




(9) 
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not contribute to the Casimir effect. The calculation is shown explicitly in 
Appendix A in 7]. 

• The above remarks were related to room temperatures. If we proceed 
to consider low temperatures, we observe that nor in this case there exist 
permittivity measurements. We thus again have to take recourse to the 
Drude relation in which, in principle, the value of v can be different from 
that above. It seems, however, that the physical significance of an altered 
value of V is only minor (a discussion on this point is given in 0). This 
is partly due to the impurities, as mentioned above, resulting in scattering 
also when the temperature is low. Important in the present context is that v 
stays finite when T — > 0, so that the behaviour is essentially as above: there 
is no contribution from the zero frequency TE mode to the Casimir force. 
Mathematically, the essential point is that 



Correspondingly, if the free energy F is drawn as a function of T at some 
fixed value of a, it turns out numerically, to a high precision, that the slope 
of the curve is zero at T = 0. This result was indicated in figure 5 in [7]; 
we intend to deal with the topic in more detail in ^Hl- That is, Nernst's 
theorem is found to be well satisfied numerically. 

• In the above argument the existence of impurities played a certain role, 
ensuring that v stays different from zero at all temperatures. Now, it is 
legitimate to ask: what about the ideal case where the metal is entirely 
free from impurities? This question, although academic, is nevertheless of 
fundamental interest. 

We may in this context recall the Bloch-Griineisen formula for the tem- 
perature dependence of the electrical resistivity p From this one may 
estimate the temperature dependence oiv = v{T^ to be 



where B = 175 K for gold. The Bloch-Griineisen argument neglects the 
effect from impurities. It is seen that z^(T) — > when T — > 0, so that this 
case becomes indeterminate. To deal with this situation, additional physical 
effects have to be drawn into consideration: 

1) One way is to include spatial dispersion, as was recently done in 
One finds by this extension of the theory practically the same results as we 





(10) 
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did above: there is only a negligible contribution to the Casimir force from 
the zero frequency TE mode. Moreover, the Nernst theorem is found to be 
satisfied, so that the presence of dissipation or a finite relaxation frequency 
in the material is not a necessity for this theorem to hold. 

2) Another approach is to take into account the anomalous skin effect 
1^ . This effect is physically due to the mean free path in the metal 
being much larger than the field penetration depth near T = 0. Again, the 
results are found to be essentially the same as above: there is no contribution 
to the Casimir force from the zero TE mode, and there is no contradiction 
with the Nernst theorem. 
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a///m 


T=l K 




T=300 K 




T=350 K 




0.16 


1144 




1127 




1124 




n 2 


508 2 




497 8 




495 7 




0.4 


38.61 




36.70 




36.35 




0.5 


16.56 




15.49 




15.30 




0.7 


4.556 




4.127 




4.052 




1.0 


1.143 




0.9852 




0.9590 




1.5 


0.2342 




0.1856 




0.1787 




2.0 


7.549 X 10- 


-2 


5.550 X 10- 


-2 


5.344 X 10" 


-2 


2.5 


3.128 X 10- 


-2 


2.176 X 10- 


-2 


2.135 X 10" 


-2 


3.0 


1.520 X 10- 


-2 


1.033 X 10- 


-2 


1.049 X 10- 


-2 


3.5 


8.252 X 10- 


-3 


5.674 X 10- 


-3 


5.990 X 10- 


-3 


4.0 


4.858 X 10- 


-3 


3.481 X 10- 


-3 


3.804 X 10- 


-3 



Table 1: The Casimir pressure between Au-Au plates versus gap width a, 
when T — {1, 300, 350} K. The pressures are given in mPa. 



a///m 


T=l K 




T=300 K 




T=350 K 




0.16 


1141 




1123 




1120 




0.2 


507.4 




496.8 




49.47 




0.4 


38.63 




36.69 




36.34 




0.5 


16.57 




15.49 




15.30 




0.7 


4.560 




4.127 




4.052 




1.0 


1.145 




0.9854 




0.9592 




1.5 


0.2345 




0.1857 




0.1787 




2.0 


7.559 X 10" 


-2 


5.551 X 10" 


-2 


5.345 X 10- 


-2 


2.5 


3.132 X 10" 


-2 


2.177 X 10- 


-2 


2.135 X 10- 


-2 


3.0 


1.522 X 10- 


-2 


1.033 X 10- 


-2 


1.049 X 10- 


-2 


3.5 


8.263 X 10- 


-3 


5.674 X 10" 


-3 


5.990 X 10" 


-3 


4.0 


4.864 X 10- 


-3 


3.481 X 10" 


-3 


3.805 X 10- 


-3 



Table 2: Same as table 1, but for Cu-Cu plates. 
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^ 1 /Tn 

CL/ jJillL 


T— 1 K 




± — OUU J-V 




J. — OtJU IV 




(J. lb 


1290 




1271 




1267 




0.2 


565.3 




553.9 




551.6 




0.4 


41.17 




39.15 




38.77 




0.5 


17.45 




16.34 




16.13 




0.7 


4.734 




4.290 




4.212 




1.0 


1.175 




1.012 




0.9853 




1.5 


0.2383 




0.1889 




0.1818 




2.0 


7.648 X 10- 


-2 


5.617 X 10" 


-2 


5.404 X 10- 


-2 


2.5 


3.160 X 10- 


-2 


2.195 X 10- 


-2 


2.150 X 10" 


-2 


3.0 


1.533 X 10- 


-2 


1.039 X 10- 


-2 


1.053 X 10" 


-2 


3.5 


8.311 X 10- 


-3 


5.693 X 10- 


-3 


6.003 X 10" 


-3 


4.0 


4.888 X 10" 


-3 


3.488 X 10" 


-3 


3.809 X 10" 


-3 



Table 3: Same as table 1, but for Al-Al plates. 



a//im 


T=l K 




T=300 K 




T=350 K 




0.16 


1143 




1125 




1122 




0.2 


507.8 




497.3 




495.2 




0.4 


38.62 




36.70 




36.34 




0.5 


16.56 




15.49 




15.30 




0.7 


4.558 




4.127 




4.052 




1.0 


1.144 




0.9853 




0.9591 




1.5 


0.2343 




0.1857 




0.1787 




2.0 


7.554 X 10" 


-2 


5.550 X 10" 


-2 


5.345 X 10" 


-2 


2.5 


3.130 X 10" 


-2 


2.177 X 10- 


-2 


2.135 X 10- 


-2 


3.0 


1.521 X 10" 


-2 


1.033 X 10- 


-2 


1.049 X 10- 


-2 


3.5 


8.258 X 10" 


-3 


5.674 X 10" 


-3 


5.990 X 10" 


-3 


4.0 


4.861 X 10" 


-3 


3.481 X 10" 


-3 


3.805 X 10" 


-3 



Table 4: Same as table 1, but for Au-Cu plates. 
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^ 1 /Tn 

CL/ jJillL 


T— 1 K 




± — OUU J-V 




J. — OtJU IV 




(J. lb 


Izlo 




line 
1195 




1191 




0.2 


535.4 




524.5 




522.3 




0.4 


39.85 




37.89 




37.52 




0.5 


16.99 




15.90 




15.70 




0.7 


4.643 




4.207 




4.130 




1.0 


1.159 




0.9986 




0.9720 




1.5 


0.2362 




0.1873 




0.1802 




2.0 


7.598 X 10- 


-2 


5.583 X 10" 


-2 


5.374 X 10- 


'2 


2.5 


3.144 X 10- 


-2 


2.185 X 10- 


-2 


2.142 X 10- 


-2 


3.0 


1.527 X 10- 


-2 


1.036 X 10- 


-2 


1.051 X 10- 


-2 


3.5 


8.281 X 10- 


-3 


5.684 X 10- 


-3 


5.996 X 10- 


-3 


4.0 


4.873 X 10- 


-3 


3.485 X 10- 


-3 


3.807 X 10- 


-3 



Table 5: Same as table 1, but for Au-Al plates. 



a//im 


T=l K 




T=300 K 




T=350 K 




0.16 


1211 




1193 




1189 




0.2 


535.0 




524.0 




521.8 




0.4 


39.86 




37.89 




37.52 




0.5 


17.00 




15.90 




15.70 




0.7 


4.646 




4.207 




4.130 




1.0 


1.159 




0.9987 




0.9720 




1.5 


0.2364 




0.1873 




0.1802 




2.0 


7.603 X 10" 


-2 


5.584 X 10- 


-2 


5.375 X 10- 


-2 


2.5 


3.146 X 10" 


-2 


2.186 X 10- 


-2 


2.143 X 10- 


-2 


3.0 


1.528 X 10" 


-2 


1.036 X 10" 


-2 


1.051 X 10" 


-2 


3.5 


8.287 X 10" 


-3 


5.684 X 10- 


-3 


5.997 X 10- 


-3 


4.0 


4.876 X 10" 


-3 


3.485 X 10- 


-3 


3.807 X 10- 


-3 



Table 6: Same as table 1, but for Cu-Al plates. 
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